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Abstract 

High resolution low frequency Raman scattering measurements from embedded AgAu nanopar- 
ticles unveil efficient scattering by harmonics of both the quadrupolar and the spherical modes. 
Comparing the experimental data with theoretical calculations that account for both the embed- 
ding medium and the resonant Raman process enables a very complete description of the observed 
multiple components in terms of harmonics of both the quadrupolar and spherical modes, with 
a dominating Raman response from the former ones. It is found that only selected harmonics 
of the quadrupolar mode contribute significantly to the Raman spectra in agreement with earlier 
theoretical predictions. 

PACS numbers: 63.22.-m,78.30.-j,78.67.-n 
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Over the past twenty years, low frequency vibrational modes from solid nanoparticles 
have emerged as original features of nanometer sized matter. These modes arise from the 
recombination of acoustic modes that turn from extended to confined upon reducing the size 
of the solid. Raman scattering is a dedicated technique to probe nanoparticle oscillation 
modes provided the used spectrometers enable to reach the low frequency range where these 
modes typically show up (a < 10 cm^ 1 ). It has proved to be a powerful tool for a quick 
and non-destructive evaluation of the sizes of the nanoparticles, especially in the case of 
metallic ones for which Raman scattering is performed in resonance with the localized sur- 
face plasmon excitations. Under these conditions, the Raman activity of selected vibration 
modes of the nanoparticles is strongly enhanced so that intense signals can be obtained 
from very diluted nanoparticle assemblies^. Metallic nanoparticles have allowed significant 
progress in the understanding of nanoparticle modes like shape anisotropy effects from Ag 
nanocolumns and nanolentils^ and elastic anisotropy effects in Au nano crystals^. Concomi- 
tantly, the theoretical description of the Raman spectra from nanoparticle vibrations has 
significantly improved over the years: while the free sphere model of H. Lamb^, combined 
with appropriate Raman selection rules^, initially provided satisfactory predictions of the 
mode frequencies^, coupling with the vibrations inside a surrounding matrix in the case of 
embedded nanoparticles have been well accounted for-^. On the other hand, quantum calcu- 
lations on free metal nanospheres, taking into account the vibration-plasmon coupling, have 
shown to be able to describe the resonant Raman spectrum in a rather complete way both in 
terms of frequencies and Raman intensities^. In particular, this latter work predicted that 
the Raman intensities from the quadrupolar mode harmonics are not a monotonic function of 
the mode index. The measurements presented in this work, which show an unusually high 
number of low-frequency Raman peaks, are an experimental proof of the non monotonic 
variation of intensity. Thanks to the observation of high order harmonics, we show that 
the interpretation of low frequency Raman data from metallic nanoparticles requires taking 
into account the coupling of the surface plasmon with the different vibrations but also the 
influence of the surrounding matrix as well as the elastic anisotropy inside the nanoparticle. 

Bimetallic AgAu nanoparticles were produced in a bulk multicomponent glass, starting 
from the incorporation of silver and gold salts to the initial oxides powder mix. After the 
colorless glass is made, prolonged annealings slightly above the glass transition tempera- 
ture generates a characteristic amber shade that arises from the surface plasmon resonance 
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FIG. 1: (color online) Absorption spectra of glasses containing Ag x Aui_ a . nanoparticles for different 
annealing times. The inset shows the SPR peak position as a function of the Ag composition, the 
dashed line is a guide to the eye. 

(SPR) of the nucleated metallic nanoparticles. The optical characterization of four samples 
annealed at 490°C over 8 H, 16 H, 32 H and 64 H shows a single surface plasmon resonance 
peak around 445 nm (Figured]). This value typically lies between that of pure silver (414 
nm) and that of pure gold (529 nm) when embedded as pure nanoparticles with comparable 
sizes in the same glass, using the same synthesis route. Upon increasing the annealing time, 
the AgAu SPR is red shifted thus reflecting a gold enrichment of the bimetallic nanopar- 
ticles. Following a classical interpretation of these datar^, one can derive the composition 
ratio x of the so-produced Ag^Aui-^ nanoparticles as a function of annealing time t a using 
a linear combination of the SPR peak positions : 

^Ag x A Ul - x (ta) = xX Ag (t a ) + (1 - x)X Au (t a ) (1) 

where \Ag(t a ) and \A u {t a ) are respectively the pure Ag and the pure Au nanoparticles 
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FIG. 2: (color online) Polarized (VV) low frequency Raman spectra (in Log-linear scale) of AgAu 
nanoparticles embedded in glasses annealed over different times. Inset: enlargement of the lower 
(E g ,T2g) doublet for the 64 H sample. 

produced in the same glass, at exactly the same annealing times. The so-derived values of 
x are found to decrease from 0.82 to 0.66 as the annealing times respectively increase from 
8 H to 64 H. 

The low frequency Raman scattering spectra were recorded with a six pass tandem Fabry- 
Perot interferometer using the 532 nm line of a continuous doubled Yttrium Aluminum Gar- 
net laser of 150 mW power. The scattered light was collected in backscattering geometry 
and the free spectral range set to 550 GHz. The finesse and optical contrast of such inter- 
ferometer are typically greater than 100 and 10 10 respectively. More details on the tandem 
Fabry-Perot interferometer can be found elsewhere^ 1 ^. 

As a general illustration, Figure [2] reports the polarized spectra (with incident and scat- 
tered light polarized both vertically, VV) of the afore characterized samples. Quite similarly 
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to recent observations^, these spectra are remarkably richer than those usually reported for 
metallic nanoparticles in the sense that they evidence more components than the usually two 
observed ones, i.e. the fundamental vibrations of both the quadrupolar mode (£ = 2, n — 1) 
and the spherical mode (I — 0, n — 1) (I denotes the angular momentum and n the harmonic 
rank). Four main bands can already be identified, some of which featuring subcomponents. 
For purposes of illustration, we will consider in the rest of the paper only the longest annealed 
sample (Ago.66Auo.34), provided that the same analysis apply equally well to all samples. 

A classical way of identifying the several components of a low frequency Raman scattering 
spectrum is to compare the polarized (W) and depolarized (VH) spectra since quadrupolar 
modes and spherical modes behave differently with regards to polarization. In a first step, 
we focus on the VH spectrum which is more simple as it only features contributions from 
the depolarized quadrupolar modes. A multi-Lorentzian decomposition of the VH spectrum 
reveals 5 components, two of them forming an intense lower frequency doublet (Figure [2] 
Inset) evidenced thanks to the high resolution Raman setup. Following recent works^, 
such doublet can be safely identified with the splitting of the fundamental mode of the 
quadrupolar vibration {£ = 2,n = 1) due to the high crystallinity of the nanoparticles (i.e. 
elastic anisotropy). Indeed, the observed frequency ratio between the lower and the larger 
components of the doublet, respectively labeled as E g and T 2g according to the irreducible 
representation of the corresponding vibrations of a face centered cubic nanocrystal, is close 
to the value expected for free Au or free Ag nano crystals^ — 1.4 us ~ 1.6 respectively). 

As suggested by previous theoretical works^ 1 ^, the higher frequency components of the 
depolarized Raman spectrum (approximately located at 150, 300 and 460 GHz) may be 
ascribed to selected high order harmonics of the quadrupolar mode. Similar predictions 
have been made recently for CdS x Sei_ x nanoparticles using different Raman scattering 
mechanisms^. In the following we compare our experimental results to Raman Scatter- 
ing Intensity (RSI) calculations based on the approach developed by Bachelier et alM. In 
this approach, the low frequency Raman spectrum of a metallic nanoparticle is generated 
from the calculation of the Raman transition probabilities using a three step scattering pro- 
cess involving the dipolar plasmon state. For the herein investigated case, this calculation 
method was extended to the case of an embedded nanoparticle, through the determina- 
tion of the continuum frequency spectrum of the nanoparticle-matrix system. A perfect 
contact between the particle and the matrix was assumed in order to apply the standard 
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boundary conditions between the two media. The full derivation of the vibration modes 
and their normalization is out of the scope of the present paper and will be described in 
details elsewhere. In addition to the excitation wavelength, the input parameters of these 
calculations are the matrix acoustic parameters (density p = 2.97 g.cm -3 ; Vl = 5020 m.s -1 
and Vt = 3010 m.s -1 , as derived from Brillouin light scattering measurements on the sam- 
ples) together with those of the bimetallic nanoparticles that were assumed to be alloys and 
thus obeying linear relationships as a function of x.— Figure [3] compares the experimental 
spectrum of the Ag .66Au .34 sample to its calculated counterpart. The overall aspect of 
the Raman spectrum is rather well reproduced in terms of band positions and relative in- 
tensities. For further discussion, the top scale of Figure [3] indicates the predictions of the 
Complex Frequency Method^ (CFM), which enables a rapid evaluation of the band posi- 
tions for embedded nanoparticles yet without yielding the Raman intensities. In order to 
identify the filiation of the Raman bands with harmonics of the quadrupolar vibration from 
a single nanoparticle, Figure [3] shows RSI spectra for the same nanoparticle but considered 
as free. In the absence of sufficiently reliable Transmission Electron Microscopy data due to 
the buried nature of the nanoparticles, the mean diameter of the nanoparticles is derived 
from the best fit of the lower frequency doublet with a single average (I — 2, n — 1) line since 
at this stage the RSI calculations are unable to account for the elastic anisotropy splitting 
(nanocrystals are assumed as isotropic nanoparticles). The derived diameter for this system 
is 23.5 nm, as confirmed from the analysis of the polarized spectrum that shows in addition 
the spherical mode contributions. 

Figure H] compares the VV and VH experimental spectra of the Ago.66Auo.34 sample, 
together with those derived from the RSI calculations. These comparisons show that the 
spherical mode contributions are rather small and they significantly overlap with those 
of the quadrupolar mode. They are better evidenced by the Raman difference spectrum 
lyy — k.Iyn, where k is a scaling factor chosen such that the (£ — 2, n — 1) peaks are 
normalized (Fig. H]left). The positions derived from the difference are found to well agree 
with the CFM predictions for the (£ = 0,n = 1,2,3) modes, using the same nanoparticle 
size as that derived from the quadrupolar mode analysis; they are equally well consistent 
with the RSI predictions of band profile changes induced by the change of polarization. 

The present observation of an usually rich low frequency Raman spectrum from embed- 
ded metallic nanoparticles, and its comparison with RSI and CFM predictions, allows to 
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FIG. 3: (color online) Experimental (symbols) and calculated (solid line) depolarized (VH) Raman 
spectra of the Ago.66Auo.34 sample with a diameter of 23.5 nm (the spectra were normalized with 
respect to the {I = 2, re = 1) peak); high frequency enlargement in the inset. The thick bars are 
the Raman intensities calculated for the same nanoparticle considered as free, labeled according to 
their n index (only harmonics with non negligible intensities are shown). Top scale arrows: CFM 
predictions for the t = 2 mode with the same embedded nanoparticle. 

capture several ingredients that are essential for a correct interpretation of similar data. 
The proper identification of the several Raman bands requires, in a first stage, to com- 
pare their behaviour as a function of polarization as they may result from the overlap of 
several contributions (Fig. H|). The detailed analysis of the depolarized spectrum, i.e. of 
the quadrupolar mode contributions, shows several important points. The first one is that 
introducing the embedding medium affects the position of the bands as displayed by the 
frequency mismatch of the observed bands (symbols in Fig. [3]) or the CFM positions (top 
arrows) with the frequencies calculated in the free case (thick bars). This effect, which is 
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FIG. 4: (color online) Left : Experimental (upper set) and calculated (lower set) of the VV (solid) 
and VH (dashed) Raman profiles of the Ago.66Auo.34 sample (intensities were normalized with 
respect to the (£ = 2,n = 1) peak). Right: Raman difference spectrum (Iyy — Iyn), experimental 
(thick line) and RSI calculated (thin line). Top scales: CFM predictions for the I = mode. 

essentially due to the relatively large mass density of the used glass, is most severe on the 
(£ = 2, n = 1) peak, which is often used for size evaluation due to its prominence. Quantita- 
tively, using a free case estimation based on the position of the (£ — 2, n — 1) band leads to 
a 19% smaller nanoparticle diameter. The second result is that from the frequency position 
matching between the experimental curve and the calculated RSI and CFM contributions, 
it comes out that the high frequency components are affiliated with selected harmonics, or 
at least harmonic groups, of the single nanoparticle vibration: the bands near 290 GHz and 
460 GHz clearly identify with respectively the n = 4, 5 and n = 8 harmonic components 
of the quadrupolar mode. Their respective Raman cross-sections are relatively consistent 
with those predicted in the free sphere case (delta functions in Fig. [3])^. The situation 
is somewhat less clearcut with the weak band at 150 Ghz which well identifies with the 
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n = 2 harmonic in terms of frequency while its intensity is considerably stronger than that 
expected from the free sphere evaluation: according to RSI calculations, it is enhanced by 
a factor 4 going from the free case to the embedded case while the other harmonics are 
enhanced by a factor 2 (its stronger intensity in the experimental spectrum is certainly due 
to its merging with the foot of the T 2g component that the RSI calculations are unable to 
account for at the present stage). This enhancement is rooted in the coupling with the ma- 
trix that significantly increases the displacement at the surface of the nanoparticle, for this 
particular harmonic, making it a potential good probe of the coupling interface. Another 
possible interpretation is that the 150 GHz band results from the redistribution of mode 
frequencies that arises from the anistropy effect, as is observed from the splitting of the 
(£ — 2,n — 1) component. These two interpretations differ from that given for a similar 
band observed from similar Au doped glasses^. 

Finally, the comparison of the RSI and CFM calculations with the experimental data 
provides a further interesting aspect related to the specificity of resonant Raman scattering 
for metallic nanoparticles. Indeed, one observes that while the frequency positions of the 
spherical mode bands agree well with those predicted by either RSI and CFM calculations 
(Fig.HJ), the CFM predictions for the quadrupolar mode harmonics are systematically slightly 
red shifted with respect to either the experimental or the RSI calculated bands (Fig. [3]). This 
is due to the dual nature of the quadrupolar modes (longitudinal vs transverse as defined 
which is also present inside a given CFM resonance. The plasmon- vibration coupling is 
therefore significant only for a subset of modes inside the frequency range defined by a CFM 
peak which can lead to a shift of the observed maxima. This situation does not occur for the 
spherical mode. As a result, the CFM and RSI peak positions match and these modes can 
be more reliably used for the size evaluation of the nanoparticles provided they are safely 
extracted from the dominating quadrupolar mode response. 

The first observation of high order quadrupolar and spherical harmonic modes from AgAu 
nanoparticles embedded in a glass has allowed us to rationalize the interpretation of low 
frequency Raman spectra from metallic nanoparticles. In accord with theoretical predictions, 
one observes that Raman scattering by the quadrupolar modes is not a monotonic decreasing 
function of the harmonic index so that only selected harmonics participate to the Raman 
spectrum with differing efficiencies. Due to stronger coupling with the plasmon excitations, 
the Raman spectrum is strongly dominated by high order harmonics of the quadrupolar 
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modes that overlap with the weaker polarized signatures from the spherical modes. In such 
case, a detailed polarization analysis of the Raman bands is therefore essential to check the 
band assignments that are uniquely based on frequency position matching 
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